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Spatial  multiplexing  in  random  wireless 

networks 

Kostas  Stamatiou,  John  G.  Proakis  and  James  R.  Zeidler 


Abstract 

We  consider  a  network  of  transmitters,  each  with  a  receiver  at  a  fixed  distance,  and  locations 
drawn  independently  according  to  a  homogeneous  Poisson  Point  Process  (PPP).  The  transmitters  and 
the  receivers  are  equipped  with  multiple  antennas.  Under  a  channel  model  that  includes  Rayleigh  fading 
and  path-loss,  and  an  outage  model  for  packet  successes,  we  examine  the  performance  of  various  spatial 
multiplexing  techniques,  namely  zero-forcing  (ZF),  ZF  with  successive  interference  cancellation  (ZF- 
SIC  or  VBLAST)  and  DBLAST.  In  each  case,  we  determine  the  number  of  streams  that  maximizes  the 
transmission  capacity,  defined  as  the  maximum  network  throughput  per  unit  area  such  that  a  constraint 
on  the  outage  probability  is  satisfied.  Numerical  results  showcase  the  benefit  of  DBLAST  over  ZF  and 
VBLAST  in  terms  of  the  transmission  capacity.  In  all  cases,  the  transmission  capacity  scales  linearly 
in  the  number  of  antennas. 


Index  Terms 

Poisson  point  process,  spatial  multiplexing,  MIMO,  outage  probability,  transmission  capacity 


1.  Introduction 

The  study  of  random  wireless  networks  has  recently  gathered  a  lot  of  momentum  in  the 
research  community,  e.g.,  see  [l]-[4].  The  main  motivation  behind  this  work  is  the  use  of 
tools  from  stochastic  geometry  in  order  to  derive  analytical  results  on  how  different  physical, 
medium-access-control  and  network  layer  parameters  affect  the  network  performance.  A  central 
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assumption  is  that  the  network  eonsists  of  a  Poisson  Point  Process  (PPP)  of  transmitters,  and 
each  transmitter  (TX)  has  a  corresponding  receiver  (RX)  at  a  given  distance.  The  justification 
for  the  widespread  use  of  this  model  is  that  it  allows  the  analytical  study  of  an  ensemble 
of  network  topologies  and  captures  the  randomness  of  the  node  locations  typical  in  networks 
without  infrastructure  such  as  ad  hoc  and  sensor  networks.  The  metric  that  quantifies  the  network 
performance  is  the  transmission  capacity,  defined  as  the  maximum  spatial  density  of  TX-RX 
links,  multiplied  by  their  rate,  such  that  a  certain  constraint  on  the  packet  success  probability  is 
satisfied  [2],  [5].  Assuming  that  the  channel  -  consisting  of  fading  and  interference  -  is  constant 
during  a  packet  slot  an  outage  model  may  be  employed  for  packet  successes,  i.e.,  a  packet  is 
successfully  received  if  the  mutual  information  of  the  channel  realization  is  greater  than  the 
desired  information  rate.  This  translates  to  a  requirement  that  the  signal-to-interference-and- 
noise-ratio  (SINR)  is  larger  than  a  predetermined  threshold. 

In  the  context  described  above,  and  provided  that  the  TX  and  the  RX  are  equipped  with 
multiple  antennas,  the  objective  of  this  paper  is  to  shed  light  on  the  performance  of  certain 
spatial  multiplexing  techniques  which  require  channel  knowledge  at  the  RX  side^;  zero-forcing 
(ZF),  ZF  with  successive  interference  cancellation  (ZF-SIC,  also  known  as  VBLAST  [6],  [7]) 
and  DBLAST  [7],  [8].  The  performance  of  these  techniques  is  well  understood  for  the  fading 
and  additive  noise  channel  [7],  but  not  so  in  a  network  environment  where  the  interfering  nodes 
are  randomly  placed. 

A.  Related  work 

The  outage  probability  and  transmission  capacity  for  different  spatial  diversity  techniques  and 
single-antenna  transmission  was  evaluated  in  [9].  One  of  the  main  results  of  this  work  was  that, 
in  the  small  outage  probability  regime,  for  maximal  ratio  combining  (MRC)  the  transmission 
capacity  scales  as  where  N  is  the  number  of  RX  antennas  and  b  is  the  path-loss  exponent. 
The  authors  in  [10]  considered  various  multiple-input  multiple-output  (MIMO)  techniques,  in¬ 
cluding  spatial  multiplexing,  as  a  component  of  a  physical  layer  that  employs  frequency  hopping 
and  coding  in  combating  interference.  They  arrived  at  similar  scaling  laws  to  [9]  regarding  the 
network  throughput  and  the  expected  progress,  albeit  from  a  different  analytical  path.  Multiple- 

'TMs  assumption  is  made  for  the  sake  of  simplicity  of  the  communication  protocol,  as  feedback  to  the  TX  is  not  required. 
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antenna  transmission  with  perfect  channel  knowledge  at  the  transmitter  was  studied  in  [11] 
and  the  optimal  number  of  spatial  modes,  in  terms  of  maximizing  the  transmission  capacity 
for  a  given  density,  was  illustrated.  More  recently,  multi-user  techniques  such  as  interference 
cancellation  and  space-division  multiple-access  have  been  considered  in  [12]-[14].  Specifically, 
in  [13],  it  was  shown  that  optimally  selecting  the  number  of  cancelled  nearby  interferers  results 
in  a  linear  scaling  of  the  transmission  capacity  with  N,  under  single- antenna  transmission. 

B.  Contributions 

We  first  consider  single-antenna  (or  single-stream)  transmission  and  revisit  the  performance 
analysis  of  MRC  in  a  Poisson  field  of  interferers,  deriving  a  compact  analytical  expression  for 
the  outage/success^  probability.  It  is  shown  that  N  RX  antennas  provide  an  approximate  gain 
of  in  terms  of  spatial  contention,  i.e.,  the  rate  of  increase  of  the  outage  probability  as  a 

function  of  the  transmitter  density,  when  the  latter  is  zero.  This  result  provides  an  alternative 
interpretation  to  the  scaling  law  derived  in  [9]. 

We  then  turn  our  attention  to  multiple-antenna  (or  multiple-stream^)  transmission  and,  em¬ 
ploying  our  findings  for  the  single-stream  scenario,  derive  exact  expressions  and  approximations 
to  the  outage  probability  for  ZF,  VBLAST  and  DBLAST.  The  optimal  number  of  streams  such 
that  the  transmission  capacity  of  the  network  is  maximized  is  determined  for  each  of  these 
techniques  in  the  small  outage  probability  regime.  The  trade-off  lies  in  the  fact  that,  introducing 
more  streams  can  potentially  boost  the  information  rate  of  each  link,  but  also  increases  the 
interference  level  in  the  network.  For  DBLAST  specifically,  it  is  shown  that,  for  6  >  4,  it  is 
optimal  to  use  all  transmit  antennas,  while,  for  6  <  4,  the  number  of  streams  must  be  judiciously 
chosen  such  that  the  optimal  trade-off  is  achieved.  Numerical  results  indicate  that  the  benefit 
of  DBLAST  over  ZF  and  VBLAST  is  significant  in  terms  of  the  transmission  capacity.  For  all 
spatial  multiplexing  techniques,  provided  that  the  number  of  streams  is  optimally  chosen,  the 
transmission  capacity  scales  linearly  in  the  number  of  antennas. 

^The  terms  “outage”  and  “success”  probability,  since  complementary,  are  used  interchangeably  throughout  the  paper. 

^In  this  paper,  the  number  of  packet  streams  is  equal  to  the  number  of  active  TX  antennas.  Each  stream  may  be  transmitted 
on  the  same  antenna,  such  as  in  ZF  or  VBLAST,  or  across  different  antennas,  as  in  DBLAST. 
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Fig.  1.  Network  model.  The  black  circles  denote  the  transmitters  and  the  green  circles  the  corresponding  receivers  at  distance 
R.  Solid/dashed  lines  denote  useful/interfering  signals. 


C.  Paper  organization  and  notation 

The  remainder  of  the  paper  is  organized  as  follows.  In  Seetion  II  we  deseribe  in  detail  our 
system  model.  Seetion  III  is  devoted  to  the  analysis  of  the  single-stream  scenario  and  Section  IV 
covers  the  extensions  to  the  multiple-stream  case.  Our  numerical  results  are  outlined  in  Section  V 
and  Section  VI  concludes  the  paper. 

We  note  the  following  regarding  the  notation:  a  zero-mean  complex  Gaussian  random  vector 
X,  with  covariance  matrix  Q  =  E[xx^]  is  denoted  as  a;  ~  CA/'(0,  Q);  the  central  chi-square 
distribution  with  parameter  1/2  and  21,  I  E  ,  degrees  of  freedom  is  denoted  as  xl/’  Ihe  I  x  I 
identity  and  zero  matrices  are  denoted  as  1^,  Oi,  respectively;  “oc”  stands  for  “proportional  to”, 
stands  for  asymptotic  equality  and  “Ri”  denotes  an  approximate  equality. 

II.  System  Model 

The  network  consists  of  an  infinite  number  of  TXs,  each  with  a  corresponding  RX  at  distance 
R,  and  locations  {xi}  that  are  drawn  independently  according  to  a  homogeneous  PPP  11  =  {xi} 
of  density  A.  Time  is  slotted  and  transmissions  take  place  concurrently  and  in  a  synchronized 
manner  during  each  slot.  Due  to  the  stationarity  of  the  homogeneous  PPP,  the  performance  of 
any  TX-RX  link,  i.e.,  “typical”  link,  may  be  studied.  The  network  model,  within  a  disc  of  finite 
radius  around  the  typical  RX,  is  depicted  in  Fig.  1. 
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The  channel  between  each  TX-RX  pair  consists  of  constant  flat  Rayleigh  fading  and  path-loss 
according  to  the  law  r~^,  with  b  >  2  (this  requirement  ensures  that  the  interference  power  is 
finite  [15]).  Additive  noise  is  disregarded,  hence  interference  from  concurrent  transmissions  is 
the  only  cause  of  errors  in  communication"^.  The  power  from  each  antenna  is  the  same  across  all 
transmitters  and,  due  to  the  absence  of  noise,  may  obtain  an  arbitrary  value,  e.g.,  unity.  Generally, 
there  is  a  different  number  of  antennas  at  the  TX  and  the  RX;  however,  for  convenience,  we 
assume  that  N  antennas  are  available  at  both  the  TX  and  the  RX^. 

Suppose  that  M  antennas  are  employed  for  transmission,  with  M  <  N.  The  received  vector 
at  the  typical  RX  can  be  written  as 

y  =  Ylx  +  w,  (1) 

where  H  is  the  N  x  M  channel  matrix  between  TX  and  RX,  with  i.i.d.  elements  [H]„m  ~ 
CJ\f{0, 1)  ;  X  CJ\f{0,  1m)  is  the  M  x  1  symbol  vector  transmitted  by  TX;  and  w  is  the 
interference  term,  modeled  as  m  ~  CA/’(0,  zl^),  where 

=  Mm  R-^  (2) 

a:iGn\{a;o} 

is  the  total  interference  power  over  a  given  slot,  per  RX  antenna;  xq  denotes  the  location  of  the 
typical  TX  and  Ri  is  the  distance  of  the  interfering  TX  at  location  Xi  from  the  typical  RX^. 
It  is  known  that  2:  is  an  a-stable  random  variable  with  stability  exponent  a  =  2/6  [1],  [4].  Its 
moment  generating  function  (mgf)  is  given  by 

$,(s)  =  E[e-^^]=e-““,  s>0,  (3) 

where  the  parameter  c  is  defined  as  c  =  A7ri?^r(l  —  a)M°‘  and  r(a;),  a;  >  0  denotes  the  gamma 
function. 

''We  select  to  study  an  interference-limited  scenario  in  order  to  focus  on  the  effect  of  cochannel  interference  on  the  performance 
of  the  employed  physical-layer  techniques.  The  analysis  can  be  generalized  to  include  thermal  noise. 

^This  assumption  is  reasonable  in  an  ad  hoc  network,  where  a  node  can  be  a  TX  or  a  RX  at  different  times. 

®Note  that,  taking  into  account  the  fading  from  an  interferer  to  a  typical  RX,  the  interference  is  generally  correlated  across 
the  RX  antennas.  Assumming  the  interference  is  uncorrelated  is  a  worst-case  scenario,  which  simplifies  the  analysis. 
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III.  Single-antenna  transmission  (M  =  1) 

Consider  the  transmission  of  a  single  stream,  i.e.,  M  =  1  and  c  =  A7ri?^r(l  —  a).  Defining  the 
desirable  information  rate  asTZ  =  log(l-|-6'),  where  9  is  an  appropriate  signal-to-interferenee-ratio 
(SIR)  threshold,  the  sueeess  probability  eorresponding  to  (1)  is  given  by  [7] 


Ps  =  P  (log  (l  +  ^)  > 

=  P(^><^).  (4) 

where  a  =  ||H|p  is  ehi-square  distributed  with  2N  degrees  of  freedom,  i.e.,  a  ~  xIn-  The 
respeetive  outage  probability  is  Po  =  I  —  Pg- 


A.  Evaluation  of  Pg 

The  evaluation  of  Pg  requires  the  knowledge  of  the  statisties  of  the  SIR  7  =  ajz.  In  the 
following  theorem,  the  eomplementary  eumulative  distribution  (eedf)  of  7  is  derived. 


Theorem  1  Let  7  =  a/ z,  where  a  ~  X2N  ^  a-stable  random  variable  with  mgf  given 
by  (3).  The  ccdf  of  y,  F^{x),  is  given  by 


N-l 


k\  ^ '  n\ 

k=l  n=k 


where 


m=l 


k 

m 


(am)n,  /c  =  1, . . .  ,n 


and  {am)n  ^  am  . . .  {am  —  n  +  1)  is  the  falling  sequential  product. 
Proof:  By  the  definition  of  Fy{x),  we  have  that 


^  +  00 


F^{x)  =  P(a  >  xz)  =  /  Fa{xy)f^{y)dy, 


where  Fa{t)  is  the  eedf  of  a,  given  by 


N-l 


f.  ,  ,  P{N,t) 


(5) 


(6) 


(7) 


(8) 
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Substituting  (8)  in  (7),  we  obtain 


^+oo 

F^{x)  =  <l>,(a;)  +  -7  / 

n=l 


'0 


From  the  Laplaee  transform  property 


IMv"  (-1)”^^ 


it  follows  that 


N-l 


n 


F^{X)  =  ^;.{X) 


n=l 


dx^' 


(9) 


(10) 


Using  identity  0.430.1,  p.24,  [16]  for  the  derivative  of  a  eomposite  function,  after  some 
algebra,  we  obtain 


d”<F^(a;) 
dx" 


=  X  e 


on 

—n^—cx^  \  '  l^k  / 


k\ 


k=l 


[CX 


(11) 


where  (3]^  is  defined  in  (6).  Substituting  (11)  in  (10)  and  regrouping  terms  results  in 


N—l  ^  n 

F,{x)  = 

n=l  k=l 


k\ 


=  e"'"*  4-  e"“ 


^  ^  (cx“)^  ^ 


k\  ^  n\ 

k=l  n=k 


(12) 


In  order  to  arrive  at  (5),  we  now  need  to  show  that  (— >  0.  Once  again,  using  the 
identity  for  the  ^th 

derivative  of  a  composite  function,  13^  can  be  written  as  the  following 
derivative  evaluated  at  x  =  1. 


^1  = 


dn(i  _x«)^ 


dx*^ 


(13) 


X  =  1 


From  (13),  the  following  iterative  relation  can  be  proved  for  n  >  2 


n 


A”=  E 

7711=1  \  ^1 


P^PkZr- 


(14) 


By  successive  application  of  (14),  we  obtain 

71  71  —  7711  — 2 - I^T-l 

J2J2'"  E  (-ir‘/3™‘(-i)”»/3™...(-i)’”‘/3;”*  (15) 


(-!)”« 


nl 


mi  =  l  7712  =  1 


mk-i=l 


where  rrik  =  n  —  ruk-i  —  ■■■  —  mi.  However,  {—lYPi  >  0,  since,  by  (6),  {—lYPi  = 
(_l)n+iQ,(Q,  _  7^  _  _  _  _l_  7^  a  =  2/6  <  1.  Therefore,  {—1Y(3^  >  0  for  /c  =  1, . . . ,  n.  ■ 


By  the  definition  of  Pg  in  (4),  we  have  that  Pg  =  F^{9)  or 


Ps 


-ce° 


+  e 


-ce^ 


N-l 


^  m 


^  k\  ^ 


n\ 


(16) 


k=l  n=k 

We  can  see  that  Pg  is  a  product  of  the  term  (the  success  probability  for  N  =  1)  and 
a  polynomial  in  of  degree  N  —  1  and  non-negative  coefficients.  Clearly,  increasing  the 
number  of  antennas  N,  increases  the  success  probability  as  more  positive  terms  are  added  to 
the  polynomial. 

In  order  to  obtain  more  insight  into  the  effect  of  iV  >  1  on  the  success  probability,  we  evaluate 
the  spatial  contention  parameter 


V 


dPg 


dX 


5 

A=0 


(17) 


defined  in  [17]  for  single- antenna  networks  as  the  slope  of  the  outage  probability  as  a  function 
of  the  density  A,  at  A  =  0.  By  its  definition,  the  larger  p  is,  the  sharper  the  increase  of  the  outage 
probability  as  A  increases.  We  have  the  following  proposition. 


Proposition  1  In  a  network  with  single-antenna  transmission  (M  =  1),  the  spatial  contention 
parameter  p  is 


T{N  -  a) 

r(iV)r(l-a)' 

Proof:  From  the  definition  of  p  and  (16),  we  have 


1/3"  I 

^  =  — -  — 

n=l 
N-l 


A  A 


cr  cr  (-l)’^(a), 
NT  NT  2-^ 


n=l 


nl 


A  ^ 


(«)n  =  1 


n=0 


n\ 


cQa  ^  I  N_l 

-  Id  2^ 


(iV-l-n)!(-l)'^-i-"(a), 


(-1)^ 

~W^  1)! 


E 


n 

N-l 

n 


(18) 
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cd^ 


-1) 


N-\ 


A  (iV-1)! 

cr  T{N  -  a) 

X  r(iv)r(i-a) 


{a  —  l)Ar_i 


(19) 

(20) 


where  (19)  stems  from  the  binomial  identity  for  falling  sequential  products  and  (20)  is  the  result 
of  the  successive  application  of  the  gamma  function  property  r(a;  +  1)  =  a;r(a;).  ■ 

For  increasing  values  of  N,  Stirling’s  approximation  yields 

~Ar-“ 


r(]v) 

However,  it  is  easy  to  verify  that  limjv^oo 

T{N  -  a) 

m) 


(l-f)^-“-^=e-“,  so 


~  N- 


(21) 


As  a  result,  for  increasing  A^,  77  ~  'kR^9°‘N~°‘ .  In  other  words,  N  provides  an  approximate  gain 
of  N~°‘  in  terms  of  spatial  contention,  or,  equivalently,  the  antenna  array  at  the  RX  effectively 
decreases  the  SIR  threshold  6*  by  a  factor  N.  Note  that  (21)  is  quite  accurate  for  relatively  small 
values  of  N,  e.g.,  for  N  =  5,  the  error  is  of  the  order  of  10%  for  6  =  4,  and  5%  for  6  =  6. 

To  conclude  the  analysis  of  the  single-stream  scenario,  the  following  proposition  provides  an 
upper  bound  on  Pg  which  is  tight  as  A  — 0. 

Proposition  2  For  single-antenna  transmission,  Pg  is  upper-bounded  as 

Ps<ew{-d>^)=Pu{X,N)  (22) 

where  the  equality  holds  for  N  =  1.  Furthermore,  for  A  — 0,  ~  Pu(A,  N). 

Proof:  For  N  =  1,  Pg  =  so  (22)  holds  as  an  equality.  For  >  1,  we  take  the  Taylor 
series  expansion  of  Pu{X,  N)  over  A  and  compare  individual  terms  with  (16).  The  reader  may 
verify  that  it  suffices  to  prove  that 


n\ 


=  Pfc, 


(23) 


n=k  \n=l 

with  /c  =  1, . . . ,  A^  —  1.  If  k  =  1,  (23)  holds  as  an  equality.  For  k  >  1,  it  holds  that 

Af-l  Af-1 

ni=l  ^2^=1 


(24) 


10 


/ _ gn 

where,  for  eonvenienee,  we  have  defined  5”  =  — -,■  Moreover,  by  (15),  we  have  that 


■i)"A 


n\ 


n  n—mi  n—m/^_2 - mi 

EE---  E  C'5r---'5r 

mfc_i=l 


(25) 


mi  =  l  7712  =  1 

where  ruk  =  n  —  ruk-i  — - mi.  Substituting  (25)  and  (24)  in  (23),  we  ean  see  that  (23)  is  a 

true  statement.  This  is  due  to  the  fact  that  the  summation  that  gives  is  over  a  subset  of  the 
terms  that  are  summed  to  give  B^. 

Finally,  by  the  definition  of  rj,  for  A  — 0,  ~  1  —  ?7A  ~  Pu{\,  N). 


As  a  result  of  Proposition  2,  (22)  can  be  used  as  an  approximation  to  Pg  in  the  small  outage 
probability  regime. 


B.  Transmission  capacity  (M  =  1) 

We  utilize  the  results  of  the  previous  subsection  in  evaluating  the  transmission  capacity  of 
the  network,  defined  as  the  maximum  network  throughput  per  unit  area,  such  that  a  constraint 
Ps  =  1  —  e  is  satisfied  [2],  [3],  i.e., 

TC,  =  A,(l  -  e)P,  (26) 

where  the  maximum  contention  density  A^  is  determined  by  the  constraint  P^  =  1  —  e.  In  the 
small  outage  probability  regime,  e.g.,  typically,  e  <  0.1,  we  can  invoke  Proposition  2  to  derive 
the  following  approximation  to  A^ 


exp 


P«(A„iV) 


/  T{N-a) 


KnR^e 


Xe 


^  1  —  € 


^  1  —  € 

^  log-(l-e)  r(7V) 

r(]V-a)' 


From  (26)  and  (27),  an  approximation  to  the  transmission  capacity  is  thus 


Ta^iog(i  +  0)- 


1)  log(l  -  e) 

7rP20« 


m) 

T{N  -a)' 


(27) 


(28) 


As  seen  by  (28)  and  (21),  the  transmission  capacity  of  a  single-stream  system  scales  as  0(A^") 


in  the  number  of  RX  antennas  N. 
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IV.  Multiple-antenna  transmission  (M  >  1) 

A.  ZF 

We  now  turn  our  attention  to  the  case  M  >  1  -  hence  c  =  XnR^Vil  —  a)M°‘.  Assume  that 
each  packet  is  transmitted  over  the  same  antenna  during  a  slot  with  a  rate  IZ  =  log(l  +  6).  If 
ZF  is  employed  at  the  RX,  the  success  probability  for  each  stream,  Pf,  is  also  given  by  (4), 
with  the  difference  that  a  is  now  chi-square  distributed  with  2{N  —  M  +  1)  degrees  of  freedom, 
as  M  —  1  degrees  of  freedom  are  sacrificed  in  order  to  cancel  out  inter-stream  interference  [7]. 
As  a  result,  invoking  (16), 


^  {ce-r  \Pl 

^  k\  ^  n\ 

k=l  n=k 


(29) 


From  Proposition  2,  we  also  have  that  Pf  <  PfX^N  —  M  +  1). 
The  transmission  capacity  is  now  defined  as 

TCf  =  Af  (1  -  e)M7^, 


(30) 


where  the  maximum  contention  density  Xf,  for  small  e,  is  determined  by  the  constraint  P„(Af ,  N— 
M  -f  1)  1  —  e  or 


A 


zf 


log(l-e)  T{N-M  +  l) 


TiP^e^  V{N  -  M  +  1-  a)M^' 

The  transmission  capacity  is  thus  given  by 

TC.  ~  log(l  +  0) - - ______ 

Due  to  (21),  for  large  values  of  —  M,  we  have  that 

''zf  1  I  1)  lc)g(l  e)  ,  ^\aT,jl-a 


TCf  ^log(l  +  0)- 


■  {N  -  M 


(31) 


(32) 


(33) 


TCf  in  (33)  can  be  analytically  optimized^  over  the  number  of  streams  M  by  allowing 
M  G  (0,  -fcx))  and  setting  the  following  derivative  to  zero 

d 


dM 


{N  -M  +  =  0. 


(34) 


’Note  that  we  can  also  optimize  over  the  SIR  threshold  6  as  in  [17]. 
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After  simple  manipulations,  we  obtain  =  (1  —  Q;)(iV  +  1).  Since  the  constraints  Mf  <  N 
and  Mf  G  must  also  be  satisfied,  the  optimal  number  of  streams  is 

Mf  =  min{r(l-a)(iV  +  l)J,iV}  (35) 

where,  with  a  slight  abuse  of  notation  [(1  —  Q;)(iV  +  1)J  denotes  the  closest  integer  number 
to  (1  —  a){N  +  1)  that  maximizes  (33).  Note  that  (1  —  Q;)(iV  +  1)  <  N  holds  if  and  only  if 
a>  l/(iV  +  l)  or6<2(iV  +  l),  which  is  valid  for  large  N  as,  typically,  6  <  6. 

Setting  M  =  (1  —  Q;)(iV  +  1)  in  (33),  we  easily  obtain  that  TCf  cx  q;"(1  —  aY~°‘{N  +  1), 
which  implies  that  TCf  =  0(iV).  The  linear  scaling  is  the  result  of  the  appropriate  choice  of 
the  number  of  streams  such  that  the  information  rate  per  MIMO  link  is  optimally  traded  off  with 
the  amount  of  interference  introduced  to  the  network.  If,  e.g.,  M  =  1  or  M  =  iV,  (32)  reveals 
that  =  0(iV“)  and  =  0(iV^““),  respectively,  i.e.,  the  scaling  is  sublinear.  This  result 
is  reminiscent  of  the  one  in  [13];  the  optimal  contention  density  scales  linearly  in  N  for  large 
N,  only  when  the  number  of  cancelled  interferers  is  a  fraction  of  N.  Interestingly,  the  optimal 
value  of  this  fraction  is  also  1  —  a. 

B.  ZF-SIC  (VBLAST) 

Suppose  that  the  RX  employs  ZF-SIC  (VBLAST),  i.e.,  it  cancels  out  each  packet  that  has 
already  been  decoded.  The  spatial  diversity  order  corresponding  to  the  “worst”  packet,  i.e., 
the  packet  that  is  decoded  first  is  N  —  M  +  1.  Given  an  outage  constraint  e  on  this  worst 
stream,  the  maximum  contention  density  is  also  given  by  (31).  Assuming  that  perfect  interference 
cancellation  takes  place,  i.e.,  there  is  no  error  propagation,  the  transmission  capacity  of  VBLAST 
is  given  by 

M 

TCf  ^  Af  log(l  +  0)  Pf  Af ,  N-m  +  l),  (36) 

m=l 

as  the  spatial  diversity  order  corresponding  to  each  stream  progressively  increases  as  more 
streams  are  subtracted  [7].  The  summation  term  in  (36)  is  the  total  throughput  of  all  transmitted 
streams,  with  corresponding  diversity  orders  —  M  +  1, . . . ,  A^,  ordered  from  the  worst  to  the 
best. 

Due  to  the  complicated  nature  of  (36),  it  is  not  possible  to  determine  analytically  the  optimum 
number  of  streams.  The  optimization  is  performed  numerically  in  Section  V. 
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C.  DBLAST 


In  Sections  IV-A  and  IV-B,  a  packet  is  transmitted  on  the  same  antenna  for  the  duration  of 
a  slot,  thereby  experiencing  the  same  fading  conditions  across  that  slot.  In  DBLAST  [7],  [18], 
a  packet  is  separated  into  segments  which  are  transmitted  across  the  antennas  and  time,  such 
that  each  segment  experiences  different  fading  conditions.  The  segments  are  then  detected  at 
the  RX  by  ZF-SIC  and,  once  a  packet  is  decoded,  its  contribution  to  the  received  signal  is 
subtracted.  It  is  known  that  DBLAST,  in  conjunction  with  appropriate  coding,  approaches*  the 
outage  performance  of  the  MIMO  Rayleigh  channel  [7],  i.e.,  for  a  total  transmission  rate  MTZ, 
the  packet  outage  probability  is  given  by 


pdh 
^  o 


P  logdet  In  +  <  MTZ  . 


(37) 


In  the  system  model  we  are  investigating,  this  probability  has  to  be  evaluated  over  the  distribu¬ 
tions  of  H  and  2:.  A  way  to  approach  this  evaluation  analytically  is  to  recall  that  may  be 
upper-bounded  as  [18] 


<  P  ^  log  (1  +  — )  <  MTZ 


=  P  (  n  ~  j  ”  <  ^  ^  I 

\m=l  / 

where  {0^}  are  independent  chi-square  random  variables  with  ~  X2(Ar-m-i-i)-  observe 
that,  for  reasonable  values  of  6  (e.g.,  6*  =  5  —  20  dB),  an  outage  roughly  occurs  when  the 
interference  power  2:  obtains  a  large  value.  In  this  case,  the  geometric  mean  of  {1  -f  is 

approximately  equal  to  the  arithmetic  mean,  thus 

/  M  j_  \  ^  \ 

P  +  ■  (39) 

\m=l  /  \  m=l  / 

(The  accuracy  of  this  approximation  is  verified  in  Section  V.)  Moreover,  ~ 


dm  \  M 


where 


Ntot  =  J2{N-m  +  l)  = 


2NM  -M^  +  M 


In  practice,  DBLAST  suffers  from  a  rate  loss  due  to  initialization. 
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As  a  result,  an  approximation  to  the  sueeess  probability  for  DBLAST,  Pf^,  can  be  obtained 
from  (16)  as  follows 


pdb  ^  ,  c(Ar(?)“  {c{M9YY  ^ 

s 


k\ 

k=l  n=k 

)db  TD  ( \  AT  \  C'r\  o  pdb 


(41) 


By  Proposition  2,  for  A  — 0,  Pf"  Pt,(A,  Atot),  so,  under  a  constraint  Pf^  =  1  —  e,  the 
optimal  contention  density  for  DBLAST  is 

log(l  -  e)  r(Atot) 


db 


'  7rP20«  r(Atot-a)M2“ 

and  the  respective  transmission  capacity  is 


TCf  ^  log(l  +  9) 


(e-l)log(l-e)  r(Atot)M 


l-2o 


7rP20«  r(Atot-«)  ■ 

From  (43)  and  (21),  for  large  values  of  Atot,  can  be  approximated  as 

-idb  ]„„/i  I  Z)\  1)  lc)s(l  ^-a/nAT  i  i\o/\/rl-o 


TCf  ^log(l  +  0)- 


7rP20o 


2-«(2A-M  +  1)“M1 


(42) 


(43) 


(44) 


As  in  Section  IV-A,  letting  M  G  (0,  +cx))  and  setting  the  derivative  of  with  respect  to 
M  equal  to  zero,  we  obtain  that  =  (1  —  a){2N  +  1).  Under  the  constraints  M  <  N  and 
M  G  ,  the  optimal  number  of  streams  for  DBLAST  is  therefore 


=  min{f(l  -a)(2A+  1)J,A} 


(45) 


where,  as  in  (35),  with  a  slight  abuse  of  notation  [(1  —  a){2N  +  1)J  denotes  the  closest  integer 
number  to  (1  —  a){2N  +  1)  that  maximizes  (44).  Note  that  (1  —  a){2N  +  1)  <  A  is  only 
possible  if  a  >  This  implies  that,  if  6  >  4  (which  is  a  typical  value  of  h  for  ground 

propagation)  transmission  with  all  antennas  maximizes  the  transmission  capacity  if  the  network 
is  operated  in  the  small  outage  probability  regime. 

We  now  investigate  how  TC^*^  scales  with  A.  Letting  M  =  in  (44)  (but  omitting  the 
operation  [-J  for  simplicity),  we  obtain  that 


f2-“a“(l-a)^-“(2A  +  l)  (46) 

TCf  oc  I 

[2-“(A  +  l)“Ai-“ 

As  a  result,  in  both  cases,  TC^^  =  0(A).  Finally,  comparing  the  optimized  transmission  capacity 
of  DBLAST  with  that  of  ZF,  we  have  that,  for  large  A  and  a  >  1/2,  TC/'^  ~  2^““TCf ,  while. 
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Fig.  2.  Success  probability  vs.  A  for  =  4  and  M  =  1  (i?  =  20  m,  6  =  4,  0  =  6  dB). 


for  a  <  1/2,  Q,c(i^_a)i-c.TCf .  The  gain  in  both  cases  is  a  direct  consequence  of  the 

robustness  of  DBLAST  with  respect  to  the  fading,  as  the  information  in  each  packet  is  coded 
and  transmitted  across  all  the  antennas  during  a  slot. 

V.  Numerical  results 

In  this  section,  we  consider  a  network  with  default  parameter  values  i?  =  20  m,  6  =  4, 
6*  =  6  dB.  In  Fig.  2,  we  plot  the  theoretical  -  eq.  (16)  -  and  simulated  success  probability  as 
a  function  of  the  PPP  density  when  M  =  1  and  iV  =  4.  The  agreement  between  theory  and 
simulation  confirms  the  validity  of  the  analysis  in  Section  III.  We  also  plot  the  upper  bound  to 
the  success  probability  given  by  (22).  As  shown  in  Proposition  2,  the  bound  becomes  tight  for 
values  of  the  success  probability  greater  than  0.8  (or,  as  the  PPP  density  becomes  progressively 
smaller). 

In  Fig.  3,  the  theoretical  and  simulated  success  probability  of  ZF  and  DBLAST  are  plotted 
vs.  the  density  of  the  PPP  for  a  system  where  M  =  3  antennas  are  employed  in  each  TX.  The 
agreement  between  theory  and  simulation  is  once  again  very  satisfactory,  which,  in  the  case  of 
DBLAST,  confirms  the  validity  of  the  approximations  in  Section  IV-C. 
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Fig.  3.  Success  probability  vs.  A  for  =  4  and  M  =  3  (R  =  20  m,  b  =  4,  9  —  6  dB). 


M 


Fig.  4.  Transmission  capacity  vs.  M  for  =  4,  8  (e  =  0.1,  i?  =  20  m,  b  =  3,  0  =  6  dB).  The  optimal  number  of  streams  for 
DBLAST  is  3  when  N  —  4,  and  6  when  N  =  8.  These  numbers  are  in  accordance  with  (45).  In  the  case  of  ZF,  the  optimal 
number  of  streams  is  1  when  N  —  4,  and  3  when  N  —  8.  These  numbers  are  also  in  accordance  with  (35).  To  avoid  cluttering 
the  figure,  DBLAST  is  denoted  as  DB  and  VBLAST  as  VB. 
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M 


Fig.  5.  Transmission  capacity  vs.  M  for  =  4,  8  (e  =  0.1,  7?  =  20  m,  &  =  4,  0  =  6  dB).  To  avoid  cluttering  the  figure, 
DBLAST  is  denoted  as  DB  and  VBLAST  as  VB. 


Fig.  4  shows  the  dependence  of  the  transmission  capacity  on  the  number  of  transmitted  streams 
for  the  three  MIMO  techniques  considered  in  Section  IV.  The  total  number  of  antennas  takes 
two  values  TV  =  4,  8,  the  propagation  exponent  is  6  =  3  and  a  constraint  e  =  0.1  is  placed  on  the 
outage  probability.  The  DBLAST  transmission  scheme  results  in  higher  transmission  capacity 
compared  to  VBLAST  or  simple  ZF.  Moreover,  the  gain  between  VBLAST  and  simple  ZF  is 
marginal,  which  is  attributed  to  the  fact  that,  with  VBLAST,  the  maximum  contention  density 
is  still  determined  by  the  subchannel  with  the  smallest  diversity  order. 

In  Fig.  5,  the  propagation  exponent  takes  the  value  6  =  4.  As  predicted  in  Section  IV-C, 
activating  all  the  TX  antennas  maximizes  the  transmission  capacity  for  DBLAST.  In  the  case  of 
ZF,  the  optimal  number  of  streams  is  dictated  by  (35).  Overall,  in  both  Fig.  4  and  Fig.  5,  the 
agreement  between  theory  and  simulation  is  satisfactory. 

In  Fig.  6  and  Fig.  7,  the  transmission  capacity  and  the  respective  -  optimal  -  number  of  streams 
are  plotted  vs.  N  for  DBLAST,  ZF  and  MRC  and  an  outage  probability  constraint  e  =  0.01. 
As  predicted  in  Section  IV,  in  the  case  of  DBLAST  and  ZF,  and  optimally  selected  M,  the 
transmission  capacity  scales  linearly  in  N,  while,  in  the  case  of  MRC,  it  scales  as  =  y/N . 
At  N  >  3,  DBLAST  provides  a  capacity  gain  of  approximately  1.4  compared  to  ZF,  which  is 
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Fig.  6.  Transmission  capacity  vs.  N  (e  —  0.01,  i?  =  20  m,  &  =  4,  0  =  6  dB).  If  all  TX  antennas  are  activated,  the  transmission 
capacity  of  ZF  is  lower  than  the  transmission  capacity  of  MRC,  even  though  the  scaling  in  both  cases  is  Setting 

M  =  N  and  6  =  4  in  (32)  we  can  also  see  that  «  rjiyh)  ’  "'hich  is  confirmed  by  the  plot. 


Fig.  7.  Optimal  number  of  streams  vs.  N  (e  =  0.01,  i?  =  20  m,  6  =  4,  0  =  6  dB). 
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in  agreement  with  the  gain  predicted  at  the  end  of  Section  IV-C.  Moreover,  it  is  observed 
that,  for  N  <  3,  MRC  and  ZF  result  in  approximately  the  same  transmission  capacity. 

VI.  Concluding  Remarks 

In  this  paper,  we  conducted  a  study  of  a  multiple-antenna  single-hop  random  network,  where 
the  locations  of  the  transmitters  are  determined  according  to  a  homogeneous  PPR  Assuming 
channel  knowledge  at  the  RX  only  and  that  interference  from  concurrent  transmissions  is  regarded 
as  noise,  we  first  evaluated  the  outage/success  probability  for  single- antenna  transmission  and 
MRC  at  the  RX.  We  then  used  the  results  and  insights  from  this  analysis  in  order  to  evaluate 
the  outage  performance  and  the  corresponding  transmission  capacity  of  MIMO  techniques  such 
as  ZF,  VBLAST  and  DBLAST.  We  determined  the  optimum  number  of  streams  such  that  the 
transmission  capacity  of  the  network  is  maximized  in  the  small  outage  probability  regime  and 
quantified  the  capacity  gain  of  DBLAST  over  ZF  and  VBLAST. 

In  conclusion,  our  results  shed  light  on  how  MIMO  techniques,  which  are  well  understood 
in  the  single-user  context,  affect  the  capacity  of  a  random  wireless  network.  At  the  heart  of  the 
analysis  and  the  resulting  design  guidelines  lies  the  PPP  geometric  model,  which  allows  us  to 
take  into  account  the  randomness  in  the  locations  of  the  interfering  TXs  in  the  statistics  of  the 
interference  power  seen  at  the  typical  RX. 
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